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ABSTRACT 

There is increasing evidence that AGN feedback is important in the energetics of cooling 
flows in galaxies and galaxy clusters. It is possible that in most cooling-flow clusters 
radiative losses from the thermal plasma are balanced, over the cluster lifetime, by 
mechanical heating from the central radio source. We investigate the implications of 
the variability of AGN mechanical luminosity L m on observations of cooling flows and 
radio galaxies in general. 

It is natural to assume that £ = ln(L m /Lx) is a Gaussian process. Then L m will be 
log-normally distributed at fixed cooling luminosity Lx, and the variance in a measure of 
L m will increase with the time-resolution of the measure. We test the consistency of these 
predictions with existing data for cooling flows and radio galaxies. These tests hinge on 
the power spectrum P(u>) of the Gaussian process £(t). General considerations suggest 
that P is a power law P ~ lj~ /3 . Long-term monitoring of Seyfert galaxies combined 
with estimates of the duty cycle of quasars imply that (5 ~ 1, which corresponds to 
flicker noise. The power spectra of microquasars have similar values of (3. We combine 
a sample of sources in cooling flows that have cavities with the assumption that the 
average mechanical luminosity of the AGN equals the cooling flow's X-ray luminosity. 
Given that the mechanical luminosities are characterized by flicker noise, we find that 
their spectral amplitudes luP(lu) lie between the estimated amplitudes of quasars and 
the measured values for the radio luminosities of microquasars. 

The model together with the observation that powerful radio galaxies lie within a 
narrow range in optical luminosity, predicts the luminosity function of radio galaxies. 
Both the shape and the normalization of the predicted function are in agreement with 
observations. Forthcoming radio surveys will test the prediction that the luminosity 
function turns over at about the smallest luminosities so far probed. 

Key words: galaxies: clusters: general - cooling flows - galaxies: jets - radio contin- 
uum: galaxies - galaxies: luminosity function, mass function 



1 INTRODUCTION 

In the cores of the majority of rich galaxy clusters the ra- 
diative cooling time of the hot, X-ray emitting intraclus- 
ter medium (ICM) is significantly shorter than the Hubble 
time. In the absence of heat sources, the gas will cool and 
flow towards the bottom of the cluster potential well: hence 
the name 'cooling-flow' clusters. In fact, recent X-ray obser- 
vations - in particular XMM-Newton high-resolution spec- 
troscopy - have ruled out the presence of significant cooling 
flows (Kaastra et al. 2001; Molendi & Pizzolato 2001; Peter- 
son et al. 2001; Tamura et al. 2001; Kaastra et al. 2004, and 
references therein): the radiating gas in the cores of these 
clusters is heated somehow. 

The most promising heating mechanism is feedback from 
the Active Galactic Nucleus (AGN) hosted by the central 
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galaxy of the cluster. The basic idea of this AGN/cooling-flow 
scenario is that the cooling ICM in the cluster core is heated 
by outbursts from the central AGN: an unsteady equilibrium 
may be achieved, in the sense that the energy input from the 
AGN balances the radiative losses of the ICM over the life- 
time of the cluster. Different models have been proposed in 
which the AGN heats the thermal plasma with either me- 
chanical (Binney & Tabor 1993; Tabor & Binney 1995) or 
radiative feedback (Ciotti & Ostriker 1997, 2001; Ostriker & 
Ciotti 2005). Although both processes might be effective to 
some extent and possibly cooperate, here we focus on me- 
chanical heating, whose relevance to many observed systems 
is apparent (e.g. Blanton et al. 2004 and Omma & Binney 
2004 and references therein). 

The cooling luminosity Lx can be determined directly 
from X-ray observations (e.g. Peres et al. 1998). By contrast, 
the mechanical power L m of the radio sources is not only 
much harder to determine observationally (e.g. Bicknell 1995; 
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Eilek 2004) , but it is expected to be time- variable (e.g. Ulrich, 
Maraschi & Urry 1997). If L m varies with time, how can it be 
compared with the cooling luminosity in observed systems? 
This paper addresses this question. 
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2 CONCEPTUAL FRAMEWORK 

The mechanical luminosity of an AGN is expected to fluc- 
tuate on time-scales that range from months to several gi- 
gayears, but from the point of view of the symbiosis between 
the AGN and its cooling flow, only time-scales in excess of 
~ 10 Myr are of interest. It is clearly impossible to probe such 
time-scales by monitoring individual AGN. Instead we must 
draw inferences by studying the demography of AGN and 
cooling flows. Such studies are still in their infancy, but it is 
nevertheless worthwhile to develop a methodology that they 
can employ, and to see what conclusions this methodology 
and currently available data lead us to. 

Cooling flows vary in size from those associated with 
large clusters of galaxies, such as Abell 1795, to those asso- 
ciated with individual elliptical galaxies, such as NGC 4472. 
We start by scaling out the variations in the mechanical lu- 
minosities L m of AGN that simply reflect variations in size 
rather than temporal variability, by focusing on the dimen- 
sionless variable 



(1) 



where Lx is the X-ray luminosity of the cooling flow, which 
is expected to vary at most slowly and weakly with time. The 
logarithm of L m , 



I = In (L m /L x ) , 



(2) 



is a continuous random variable defined on the interval 
(—00,00). Throughout physics one models such a variable 
through the Fourier representation 



*(«) = / 



(3) 



The complex numbers are random variables. The sim- 
plest conjecture is that these variables are statistically inde- 
pendent of one another, and have phases that are uniformly 
distributed in (0, 27r). In this case £(t) is a Gaussian ran- 
dom process, and it is completely characterized by the power 
spectrum P{y})= {\A^\ 2 ). The probability density of the lu- 
minosity at any given time t is given by 



P(L m )dL„ 



271-er 



-exp 



2<7 2 



At, 



(4) 



where a 2 — ((£ — {£)) 2 ) is the variance of i. It is easy to show 
that with this model the expectation value of the nth power 



of L n 



(O = exp ( n{£) +n 2 y 



(5) 



The appearance of a 2 in the exponential signals the skew- 
ness of the probability distribution Q. The larger the value 
of n, the more important this skewness is in the sense that 
rare excursions to high L m more strongly dominate (Li^). A 
measure of the importance of skewness is the probability that 
an individual measurement of L m will return a value smaller 
than (L m ). One can readily show that 



P(L m < (L m » = ±[1 + erf (a/2 3 / 2 )]. 



(6) 
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Figure 1. The probability for a log- normal distribution that a 
single measurement yields a value smaller than the expectation 
value. 



Figure^plots this function of a. When a > 3, there is a high 
probability of finding the luminosity to be smaller than its 
mean value because much of the total energy output comes 
during infrequent outbursts. 

We are interested in the properties of the light-curve 
observed with a finite time-resolution, i.e., averaged over a 
given time At. We define this 'coarse-grained' light-curve as 



p 00 

*At(t)= / dt' FAt(\t' -t\)£(t'), 
Jo 



(7) 



where the filter function FAt{t) satisfies J dtFAt = 1 and is 
non-zero in a region of width At around the origin. ^At(t) 
has average 



(£ A t) = lim - / dt'lAt(t') = {£), 



t — >oo t 

and variance 1 

a 2 At = lim - / di(i A t(t')- {£)f 

t^oo t Jq 



dto Fi t (uj)P(uj) 



(8) 



(9) 



where i<At(w) is the Fourier transform of the filter FAt{t). 
Approximating the window function F At (tu) with the Heavi- 
side step function 9(1 — w/wAt), where u>At = 2n/At, we get 



2 

&At 



du>P(u 



(10) 



In words, contributions to the variance on time-scales shorter 
than At vanish in the averaging process. Thus, for any power 
spectrum P(w), a At i s a non-increasing function of At, with 
limAt^oo <y 2 At ~ and limAt^o &At ~ a2 '- the variance is 
smaller when averaging over longer At. 



2.1 The power spectrum 

An AGN involves an enormous dynamical range of time- 
scales, from the few days characteristic of the dynamical time 
at the event horizon to the hundreds of megayears required 
for the development of a cooling catastrophe in the core of 



1 Note the strict similarity between the variance of the coarse- 
grained light-curve cr\ t and the mass variance a 2 M usually adopted 
in cosmology to describe mass density fluctuations on mass scales 
larger than M. 
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the cooling flow. In view of this great span of time-scales, 
one might expect the power spectrum of the AGN to be a 
power law in frequency. So we assume that £(t) has a power 
spectrum in the form 



P (lo/lo ) 13 ifcun 



< id < Wn 







otherwise, 



(11) 



where Po is the power at some reference frequency too , iJmin = 
27r/t m ax and a; max = 2ir/t m in- t m in and i max are, respectively, 
the maximum and minimum time-scale of variability. 

There is no difficulty in principle in deriving equations 
for a general power-law index /3, but we shall argue that j3 
must lie close to unity, and the discussion is greatly simplified 
if we specialize to j3 = 1 at this point. The case /3 = 1 
is associated with 'flicker noise', which empirically arises in 
diverse physical phenomena that range from electrical noise, 
to wind speeds, to earthquake magnitudes (e.g. Press 1978; 
Milotti 2002) . In case of flicker noise the relation between the 
variance and the power spectrum is 



2 T-, , ^Ai n , imax 

"At = Po^o m = P LO O m — r— 

(Jmin At 



(12) 



The amplitude of flicker noise is conveniently characterized 
by the dimensionless quantity uP(lo), which is independent 
of LO. 



2.1.1 Calibration with quasars 

The variability of X-ray emission from AGN has been studied 
on time-scales shorter than ~ 10 yr. When the X-ray power 
spectra 2 of Seyfert galaxies at frequencies 10 -8 — 10 -5 Hz are 
fitted with a power-law, one finds /3 ~ 1 (for instance, (5 ~ 
1.05t°;2o for the Seyfert galaxy NGC4051), with amplitudes 
P(lu) ~ 0.01 - 0.04 (Uttley, McHardy & Papadakis 2002; 
McHardy et al. 2004). Hence, the only power-law indices that 
we can consistently adopt are (5 ~ 1. 

There is a vast difference between the longest time-scales 
on which the power spectra have been empirically deter- 
mined, and the time-scales > 10 7 yr of interest here. Hence 
to extrapolate blindly from the observed < 10 yr fluctuations 
out to 10 7 yr ones is a very hazardous proceeding. Fortu- 
nately, independent estimates of the duty cycle of quasars 
provides a degree of reassurance. Ciotti & Ostriker (2001) 
define the duty cycle of a source with light-curve L(t) as 



{I?) 



(13) 



where the second equality follows from equation Var- 
ious arguments indicate that quasars, over the lifetime of 



their host galaxies have duty cycles 5 ~ 10 



10 , with 



episodic lifetime (i.e., duration of individual major outbursts) 
> 10 4 yr (Haehnelt, Natarajan & Rees 1998; Yu & Tremaine 
2002; Martini & Schneider 2003; Haiman, Ciotti & Ostriker 
2004). By equation 1131 . these duty cycles require vari- 
ances in the range ~ 4.6 — 6.9. Inserting these variances 
and the time-scale range 10 4 — 10 10 yr into equation (1121 . 
we find that the associated flicker spectra have amplitudes 
Poujo ~ 0.3 — 0.5. These values are larger by a factor ~ 30 

2 The authors consider the power spectrum of L(t) normalized to 
(L) 2 , while here we refer to the power spectrum of i(t) = lnL(t): 
however, the two have the same slope and amplitude as long as 
the variance of L(t) is small with respect to (L) 2 . 



than the amplitudes measured by Uttley et al. (2002) for 
Seyfert galaxies at frequencies higher by >10 6 . If f5 dif- 
fered significantly from unity, we would expect the ampli- 
tudes to be many orders of magnitude different, so we must 
have 13 ~ 1. 



2.1.2 Calibration from microquasars 

Microquasars are radio-jet emitting X-ray binaries made up 
by a stellar-mass black hole (BH) and a companion star, and 
seem in many respects small-scale versions of AGN (Mirabel 
& Rodriguez 1999). It is widely accepted that on small to 
intermediate scales these objects are identical to AGN scaled 
down in both length- and time-scales by the ratio of the 
black holes' Schwarzschild radii r g = 2GMbh/c 2 . Hence, we 
can study the temporal variability of microquasars on hu- 
man time-scales to get information about the variability of 
AGN on time-scales much longer than those that can be di- 
rectly probed. It seems likely that the scaling holds at out 
to time-scales 3 ~ 10 5 yr, which for an AGN black-hole mass 
10 7 — lO 8 A/0 corresponds to a readily observed frequency 
10 _5 Hz in microquasars. Thus from studies of microquasars 
we can estimate the amplitude of the AGN power spectrum 
at frequencies corresponding to 10 yr, and then use our 
knowledge that /3 cannot significantly differ from unity to 
extrapolate it to frequencies lower by a factor ~ 100. 

Most studies of the low-frequency power spectrum of mi- 
croquasars focus on X-ray emission. Reig, Papadakis & Ky- 
lafis (2002) find that the power spectrum of the microquasar 
GRS 1915+105 is well fitted by a power law f equation 1111 
with index /3 ~ 1 in the frequency range 10 -7 to 10 -5 Hz. 
A power spectrum with this slope, over the same frequency 
range, has been also found for the X-ray light curve of Cygnus 
X-3 by Choudhury et al. (2004). 

However, we are interested in the variability of the ki- 
netic power of the jet, which is expected to be traced by radio 
emission rather than by X-rays. The radio variability of mi- 
croquasars has not been investigated as much as their X-ray 
variability. A spectral analysis was carried out by Fiedler at 
al. (1987) for the radio (2.25 and 8.3 GHz) light curve of 
the peculiar microquasar SS 433: they found that its power 
spectrum is consistent with a power law with index /3 ~ 1 
at frequencies 10 -7 to 10 -6 Hz (see also Revnivtsev et al. 
2005). Nipoti, Blundell & Binney (2005) show that at these 
low frequencies the slopes of the radio power spectra of the 
microquasars Cygnus X-3 and GRS 1915+105 are similar to 
those found in X-rays. 

The X-ray light-curve of GRS 1915+105 has power spec- 
trum with amplitude Pqloq ~ 0.02 at frequencies 10 -7 to 
10 Hz (Reig et al. 2002). An amplitude smaller by a factor 
of ~ 10 has been measured for the X-ray power spectrum 
of Cygnus X-3 in the same frequency range (Choudhury et 
al. 2004). However, there is evidence that the amplitude of 
the variability of the radio emission of these powerful micro- 
quasars is larger than that of their X-ray emission (Nipoti 
et al. 2005). The radio emission of Cygnus X-3 (at 2.25 and 



3 For the accretion disk of a ~ 10 s Mq BH, the characteristic 
viscous time-scale at ~ 10 3 r , g is roughly 10 5 yr. Longer time-scales 
correspond to larger radii, at which AGN accretion disks might be 
truncated by self-gravity (e.g. Burderi, King & Szuszkiewicz 1998; 
Goodman 2003). 
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Figure 2. Possible realization of the radio light-curve (normalized 
to the averaged radio luminosity) of an AGN. The light-curve is 
obtained from the daily-averaged 2.25 GHz radio light-curve of 
Cygnus X-3 over two decades by rescaling in time by a factor 10 8 . 



8.3 GHz) has been monitored by the GBI-NASA monitor- 
ing program (Waltman et al. 1995, and references therein). 
We find that the daily-averaged 2.25-GHz radio light-curve 
of Cygnus X-3 has logarithmic variance a 2 ~ 0.96 over the 
period 1982-2001. If we adopt t max = 19 yr, At = 1 day for 
the boundaries of the frequency range sampled, equation 1121 1 
implies that the power spectrum of the radio emission has 
amplitude Poi^o — 0.11, about two orders of magnitude 
larger than that found in X-rays. Similarly, we estimated the 
amplitude of the radio power spectrum of the microquasar 
GRS 1915+105, using its 2.25-GHz light-curve obtained by 
the GBI monitoring program over about 6yr (see Foster 
et al. 1996). We found that GRS 1915+105 has amplitude 
Pqluq ~ 0.12 at 2.25 GHz, also much larger than in X-rays. 
Thus, observations of microquasars suggest Poojq ~ 0.1 as 
a plausible value of the amplitude of the very low-frequency 
power spectrum of radio sources. Given the substantial uncer- 
tainies involved, this estimate is in remarkably good agree- 
ment with the estimate Poloo ~ 0.4 that we obtained by 
considering the duty cycle of quasars. 

These considerations suggest that we can obtain an im- 
pression of the long-term variability of AGN by rescaling the 
radio-frequency light curve of a microquasar. Fig. [5] is a plot 
of the daily- averaged 2.25 GHz light-curve of Cygnus X-3 
rescaled in time by 10 8 . The light curve is characterized by 
short powerful outbursts, alternating with periods of quies- 
cence. 



3 APPLICATION TO COOLING FLOWS 
3.1 Systems with cavities 

In general, it is not feasible to measure directly the mechani- 
cal power of radio sources. An exception is M87 - the central 
galaxy of the Virgo cluster - for which high-resolution obser- 



vations of the jet and the very inner hot atmosphere are avail- 
able (Owen, Eilek & Kassim 2000; Di Matteo et al. 2003). In 
the case of the central radio sources in galaxy clusters L m 
can be estimated from depressions in the X-ray emission of 
the surrounding ICM. These 'cavities' are interpreted as bub- 
bles produced in the interaction between the jets of the radio 
source and the thermal plasma. The energy transferred to the 
ICM by a bubble expanding adiabatically is (e.g. Churazov 
et al. 2000, 2002) 



E h = 



7 



7- 



- P V, 



(14) 



where V is the volume of the bubble, p is the pressure of 
the surrounding ICM, and 7 is the adiabatic index of the 
gas m the bubble. Thus, E h = \pV for a non-relativistic 
gas (7 = |), and E\> — 4pV for a relativistic fluid (7 = |). 
Given a characteristic time Atb associated with the bubble, 
an estimate of the average kinetic luminosity of the AGN 
over Atb is then 



E b 
At b ' 



(15) 



The definition of Atb is controversial: if we consider Atb 
as the age of the bubble, this is of order the time t Ca required 
for the cavity to rise at the speed of sound from the centre to 
its current location (Gull & Northover 1973; Churazov et al. 
2001). Alternatively, the age of the bubble can be defined as 
the time required for the ICM to refill the displaced volume 
(McNamara et al. 2000). However, Atb is not necessarily the 
age of the bubble: it could be the recurrence time t ICC of 
the radio source, defined as the time elapsing between two 
subsequent outbursts in which the jet blows bubbles. As we 
are looking for the characteristic time to be associated with 
this measure of L m (i.e., the time over which we average in 
measuring from the bubbles), it is natural to identify Atb 
with t roc , so in the rest of the paper we assume Atb = t Icc . 

Birzan et al. (2004, hereafter B04) present radio and 
X-ray data for a sample of 18 systems with cooling flows 
that show evidence of cavities in their diffuse X-ray emission. 
The sample consists of 16 galaxy clusters, the galaxy group 
HCG 62, and the elliptical galaxy M84. The authors provide 
for each of these objects the cooling luminosity (Lx), the 
monochromatic radio luminosity at 1.4 GHz (L1.4), the inte- 
grated radio luminosity between 10 MHz and 5 GHz (I/ ra d)i 
and measures of pV for each cavity. This sample is suitable for 
comparison with the AGN/cooling flow model, because for 
each source different measurements of the mechanical power 
(averaged on different time-scales) are available. 

From measures of pV we derive L m and L m = L m /Lx us- 
ing equation 1151 . given the time-scale Atb = ircc- In general, 
we do not know the recurrence time for individual sources. 
In the Perseus cluster, in which two generations of bubbles 
are observed (Fabian et al. 2000) , the recurrence time-scale is 
of the order of 100 Myr. Estimating t ICC for other sources, in 
which just one generation of cavities is observed, is problem- 
atic, but statistical arguments suggest t lac ~ few x 100 Myr 
as a typical value (B04) . So we assume that Atb is the same 
for all sources, and we parameterize L m as 



/b 



pV 
10 8 yr ' 



where 



7 ( At 



7-1 V 10 8 yr 



(16) 



(17) 
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Figure 3. 90 and 95 per cent confidence levels in the space of the 
parameters fb = [7/(7 — l)]/(Att,/10 8 yr) and variance a 2 for the 
distribution of L m (measured from bubbles) for B04's sample of 18 
objects with cavities. The null hypothesis is that the distribution 
is log-normal with (L m ) = 1. 

We then use the Kolmogorov-Smirnov (KS) test to esti- 
mate how well the measured values of L m fit the predicted 
log- normal distribution @ with (L m ) = 1. Fig. |3] plots 
the 90 and 95 per cent confidence levels, derived from the 
KS test, in the space of the parameters fb and a 2 . For 
instance, if we fix 7/(7 — 1) = 4, and allow Atb in the 
range 1 < (Afb/10 8 yr) < 4, the factor fb lies in the interval 
1 < fb < 4. Such a range in fb corresponds to values of the 
variance 0-2 <cr 2 <3.9 that give distributions consistent with 
the data at the 95 per cent confidence level. 

3.2 Calibrating the radio luminosity 

The number of systems in which cavities have been accu- 
rately observed is small, while the radio luminosities of thou- 
sands of radio galaxies have been measured. It is clear that 
the radio power represents a very small fraction of the ki- 
netic power emitted by the AGN (e.g. Eilek 2004), but the 
radio luminosity may be a useful indicator of the mechanical 
power. Hence, much larger samples could be assembled if it 
were possible to establish a relation between the mechanical 
and radio luminosities of sources. 

Cavities in the X-ray emitting ICM are classified as 
radio-filled if non-thermal radio emission at a given frequency 
coincides with the X-ray-depressed area; otherwise they are 
called 'ghost' cavities 4 . If reacceleration of the synchrotron 
electrons is negligible, the age of the radio-filled bubbles will 
be of the order of the synchrotron lifetime, or smaller, while 
the age of the ghosts will be longer. Thus, from radio-filled 
cavities we can independently measure L m ~ Eb/t Cs and L1.4 
and thus determine the relation that holds between them. 

The top panel of Fig. [I] plots data for the nine objects 5 

4 This classification depends on the frequency of the radio emis- 
sion. Cavities classified as ghosts on the basis of ~ GHz radio data 
might appear as radio-filled at lower frequencies, depending on the 
break frequency of the radio spectrum. This has been observed in 
Perseus by Fabian et al. (2002). 

5 We do not include Abell 478, whose cavities are considered 
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Figure 4. The ratio pV/t Cs as a function of the monochromatic 
1.4 GHz luminosity (top) and the bolometric radio luminosity 
(bottom) for radio-filled cavities (data from B04). The solid lines 
are the best-fitting linear relations: equations 1181 and (1191 with 
M.4 = 897/(7 — 1), and K ra( j = 11-17/(7 — 1), respectively, assum- 
ing L m = E h /t Cs . 

with radio-filled cavities presented by B04. The values of pV 
are from Chandra, and t CB is from B04. L1.4 — Pi. 4 x 1.4 GHz, 
where Pi. 4 is the monochromatic 1.4 GHz power measured by 
the VLA. The straight line is the relation 

Lm = ^—7 — = K1.4I/I.4, (18) 

7 - 1 *c. 

with the best-fitting dimensionless conversion factor kia — 
897/(7-1). 

The lower panel of Fig. 0] shows the result of repeating 
this exercise using the bolometric radio luminosity L ra( j in- 
stead of the monochromatic luminosity. L ra< j is estimated by 
integrated between 10 MHz and 5 GHz, assuming constant 
spectral index for each source (B04). The straight line is the 
relation 

L m — ^rad-^rad, (19) 

with best-fitting K ra d — 11.17/(7— !)■ B04 found that a non- 
linear relation - the mechanical luminosity scaling as L^j 1 
- provides a better fit to their data (including Abell 478) 

radio-filled by B04, but are more likely to be already fading to 
ghosts, as suggested by Sun et al. (2003). Also B04 point out the 
anomalous high ratio of mechanical to radio power of this source. 
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Figure 5. Top: 90 and 95 per cent confidence levels, in the space 
of the parameters rei.4 (equation 1181 and variance a 2 4 for the 
distribution of L m (measured from L14) for B04's sample. Bot- 
tom: the same as in top panel, but in the space re r ad _cr rad' ^ or * ne 
distribution of L m measured from -L ra d- The vertical dotted lines 
correspond to the best— fitting values K1.4 ~ 356 and K ra d ~ 44 
found for the radio-filled cavities, assuming 7/(7 — 1) = 4. The null 
hypothesis is that the distributions are log-normal with (L m ) = 1. 



than the simple proportionality L m oc L ra d considered here. 
However, the sample is very small, and one can show that, 
excluding Abell 478 and the peculiar radio source Cygnus A - 
the only FR II (Fanaroff & Riley 1974) source in the sample 
- the slope of the best-fitting power-law is close to 1. We 
think that a larger sample is needed to verify the existence 
of a non-linear relation between kinetic and radio power. 

We next investigate the constraints that can be placed on 
Kia and Krad by assuming that (L m ) = Lx- We use the whole 
sample presented by B04, without any distinction between 
radio-filled and ghost cavities, and assume that the conver- 
sion coefficients kia and K ra( j are the same for all sources. We 
use the KS test to evaluate the consistency between the val- 
ues of L m obtained using trial values of these coefficients and 
the log-normal probability distribution Q with (L m ) = 1- 
Fig. shows contours of 90 and 95 percent confidence lev- 
els in the planes of a 2 versus log(/ti.4) (top) or log(K ra d) 
(bottom). The dotted line in each panel shows the value of 
k determined above from Fig. 2] under the assumption that 



7/(7 — 1) = 4. The data are consistent for a wide range of 
parameters. Larger variances a 2 are required for smaller co- 
efficients k because the current values of L m are proportional 
to k, and if these are small the assumption (L m ) = 1 requires 
that significant energy is released during excursions to high 
luminosity. For the values of k that we derived from radio- 
filled bubbles, of ad ~ 3 and erf 4 ~ 4 are favoured variances. 
These values are intermediate between the values w 1 that 
we inferred from microquasars and ~ 5 that we estimated 
from quasars. 

3.3 Dependence of a 2 on time-scale 

As discussed in Section 2, the variance a 2 associated with any 
measurement of L m should decrease with increasing timespan 
over which that measurement effectively averages L m . Using 
radio data we estimate the kinetic power averaged over the 
lifetime of synchrotron-emitting relativistic electrons, which 
depends on the radio frequency. A reference time-scale for 
synchrotron radiation at frequency v is the synchrotron cool- 
ing time 

/ B \-3/2 1/2 

^ 5 °U^g) (gs) Myr ' (20) 

where B is the modulus of the magnetic field (e.g. Carilli 
et al. 1991). t syac depends strongly on B, which is usually 
poorly constrained in observed radio sources, and for fixed 
strength of the magnetic field decreases by about one order 
of magnitude from v = 70 MHz to v — 5 GHz. Actually, t sync 
represents an upper limit on the lifetime of the radio-emitting 
electrons because adiabatic expansion losses are expected to 
be dominant over radiative losses (e.g. Blundell & Rawlings 
2000). As the effects of adiabatic expansion losses are diffi- 
cult to quantify in general, in our treatment we simply as- 
sume that the averaging time-scale is < i syn c (y), when L m is 
measured from radio luminosity at frequency v. 

Fig. plots the three ranges of acceptable values of 
a 2 that we have derived above against the characteristic 
averaging time At associated with each measurement. For 
the mechanical luminosity measured from the bubbles, we 
have 0.2<a^<3.9, and we adopt 100 < At h < 400 Myr as 
the range for the recurrence time (see Section 3.1). For the 
mechanical luminosity derived from 1.4 GHz radio data, 
1.8 < erf 4 < 6.9 and 3< Ati.4 < 36 Myr, where the range in 
Ati. 4 = i sync (1.4 GHz) has been computed from equa- 
tion H20H for B in the range (1 — 5) x 10~ 5 G - consistent 
with minimum-energy magnetic fields measured in the ra- 
dio lobes of 3C 218 (Hydra A; Taylor et. al 1990), 3C 84 
(Perseus; Fabian et al. 2002), and 3C 405 (Cygnus A; Car- 
illi et al. 1991). Considering the bolometric radio luminosity, 
we get more stringent limits on the variance, 1.9 < a 2 ^ < 4.6, 
while the corresponding Ai ra d spans the same range as Ati.4, 
because in this case the bolometric luminosity is measured 
from 1.4 GHz radiation. We assume that the relevant limits 
on the variance are those found for a 2 ^ d , because the spread 
in the slopes of the radio spectra is likely to contribute to 

„2 

CT4- 

The dashed curves in Fig. show the dependence of 
a 2 upon At that is expected for flicker noise spectra (equa- 
tion I10II with two different amplitudes, Poi^o = 0.24 and 
0.84. The model and the data are consistent for amplitudes 
that lie within this range, which overlaps the interval be- 
tween the amplitudes of the radio spectra of microquasars 
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Figure 6. Variance of the logarithm of the mechanical power of 
a radio source, as a function of At, the time over which the me- 
chanical power is averaged. The two solid boxes define the loci 
of the diagram with pairs (At^o^) and (At ra[ j,o"^ ad ) derived for 
B04's sample from observations of X-ray cavities and bolometric 
radio luminosity, respectively. The dashed curves are the upper 
and lower limit of the family of functions <j\t predicted for a 1// 
power spectrum, which are consistent with the data. The dotted 
box refers to pairs (Ati.4,(rJ 4 ) derived from monochromatic 1.4 
GHz data for the same sample. 



and those inferred from the duty cycle of quasars. Note that 
for Atr a d we used the synchrotron cooling time, which is an 
upper limit. A shorter At ra d would result in lower value of 
Pouio: for instance, reducing At ra d by a factor of 10 we get 
0.19 <P uo< 0.59. 



4 APPLICATION TO THE RADIO-GALAXY 
LUMINOSITY FUNCTION 

In the previous section we have shown that, under plausi- 
ble assumptions, observed data of central radio galaxies in 
cooling-flow clusters are consistent with the hypothesis that 
AGN heating balances radiative cooling. Here we discuss the 
consequences of extending the considered picture to cooling- 
flow systems in general. In particular, applying our model to 
luminous early-type galaxies leads to predictions on the local 
radio luminosity function. 

Radio galaxies, though spanning several orders of mag- 
nitude in radio power, represent a very homogeneous fam- 
ily of objects when the properties of their host stellar sys- 
tems are considered, being early-type galaxies with optical 
luminosity larger than L* . For instance, Magliocchetti et al. 
(2002) derived the local radio luminosity function for AGN 
and early-type galaxies of the FIRST-2dF Galaxy Redshift 
Survey (2dFGRS), finding that the optical counterparts have 
B-band luminosity characterized by a very narrow distribu- 
tion, with average A4b — 5 log h = —19.95 and scatter 0.41 
mag [h = i/ /(100 kms^ Mpc" 1 ), where H is the Hubble 
constant] . 



Figure 7. Local radio luminosity function of early-type galax- 
ies. Points with error bars are data from the FIRST-2dFGRS 
(Magliocchetti et al. 2002). The dashed line is the best-fitting 
log- normal distribution ( equation 1231 . 



Let us assume that all bright elliptical galaxies have a 
cooling core, a central supermassive BH, and that, over the 
galaxy lifetime, the radiative losses from their hot plasma are 
balanced by the mechanical energy emitted in radio outbursts 
from the central AGN ((L m ) = Lx)- The radio luminosity 
of the AGN will be strongly variable, and poorly correlated 
with the cooling luminosity of the host cooling core. However, 
under the assumption that the monochromatic 1.4 GHz radio 
luminosity traces the kinetic luminosity (L m = K1.4L1.4), for 
a sample of galaxies with the same cooling luminosity Lx , we 
would expect (L1.4) = Lx/ki.4- If ^>(Lx) is the distribution 
of the cooling luminosities of bright early-type galaxies, the 
radio luminosity function of early-type galaxies will be given 
by 



$(Li 



(21) 



where V(L\a\ Lx) is a log-normal distribution with {L1.4) = 
Lx/k-ia and variance CX1.4. The luminosity function vE'(Lx) 
has yet to be determined observationally, because of the lack 
of complete X-ray samples of early- type galaxies. As well 
known, the diffuse X-ray luminosity of early-type galaxies 
with similar optical luminosity is characterized by signifi- 
cant scatter (Fabbiano, Gioia & Trinchieri 1989; O'Sullivan, 
Forbes & Ponman 2001). Given that the optical luminosity 
is roughly the same for all sources, it is natural to model the 
distribution of Lx as 



*(Lx)dL x 



No 



lixox 



■ exp 



(^x-(^)) 2 



2ol 



d£x 



(22) 



where £x = ln(Lx), o\ i s the logarithmic variance and No is 
the number density. Substituting >&(Lx) and V(L\a\Lx) in 
equation 12 1 > we get 



$(Li. 4 )dLi. 4 = 



No 



27TCT, 



exp 



InLt. 



2ol 



dii.4, (23) 
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where £1.4 = ln(Li.4), <x* = a\ A + <j x , and mLi 4 = 
(^x) — In K1.4 — af 4 /2. Thus, our model predicts that the ra- 
dio luminosity function of early-type galaxies is log-normal 
with variance 171.4 + °x> where a\ A is the variance due to 
radio-variability of each source and a x is the variance in the 
distribution of X-ray emission of bright ellipticals. 

We test this prediction by fitting the local radio lu- 
minosity function of early-type galaxies of the FIRST- 
2dFGRS (Magliocchetti et al. 2002) with equation QZQ, leav- 
ing ctJ, L\ 4 and No as free parameters. We found best- 
fitting parameters (with 1 a uncertainties) a'l = 7.3 ± 1.0, 
LI .4 = 10 38 - 5±0 - 2 erg s" 1 , and iV = (1.2 ± 0.3) x 
10 -3 /i 3 Mpc~ 3 . Fig. 7 plots the data and error bars as re- 
ported by Magliocchetti et al. (2002) together with our best- 
fitting log-normal function. We note that the best-fit gives 
Cx.4+0x — 7-3, consistent with values of a\ A found analysing 
B04's sample of central radio galaxies in cooling cores (Sec- 
tion 3.2) plus a contribution from intrinsic scatter in X-ray 
luminosity. 

We based our analysis of the radio luminosity function 
on the assumption that most bright elliptical and SO galax- 
ies have AGN-powered radio emission. We can check this 
assumption a posteriori by comparing the best-fitting num- 
ber density TVo with the number density of luminous early- 
type galaxies. We assume that the B-band optical luminosity 
function of early-type galaxies <^>(Lb) is a Schechter function 
with slope a, characteristic luminosity L B , and normalization 
4>* . Then the number density of objects brighter than L B is 
N B (> L B )=j™dL' B 4>(L B ) = 4>*T(a+l, Lb/L b ). If we con- 
sider the best-fitting B-band luminosity function determined 
by Madgwick et al. (2002) for the early-type galaxies of the 
2dFGRS [a = -0.54T0.02, <j>* = (9.9T0.5) x 10 _3 /i 3 Mpc" 3 ], 
and we assume that all early-type galaxies with L B > L B 
are radio galaxies, we find a number density N B (> L B ) = 
(2.7 ± 0.1) x 10" 3 /i 3 Mpc- 3 . Thus, iV B (> L%) is a factor of 
~ 2 larger than TVo found from our best-fitting log-normal 
distribution. The discrepancy between the two number den- 
sities is easily explained by our arbitrary choice of L B as 
the minimum B-band luminosity of a radio galaxy, and by 
uncertainties on the faint-end of radio luminosity function. 

We conclude that available radio data are consistent 
with the proposed scenario, in which each bright elliptical 
oscillates, over its lifetime, along the radio luminosity func- 
tion, spanning several orders of magnitude in radio power. 
Interestingly, results of high-redshift radio surveys also sug- 
gest that most massive ellipticals have experienced one or 
more episodes of powerful radio activity during their lifetime 
(Rawlings 2003) . Finally, we stress that the proposed picture 
naturally explains the very large scatter in radio luminosity 
for radio galaxies with the same BH mass (Lacy et al. 2001; 
McLure et al. 2004, and references therein) as a consequence 
of the strong time- variability of their radio emission (see also 
Nipoti et al. 2005). 



5 SUMMARY AND CONCLUSIONS 

Since the cosmic density of matter is very inhomogeneous 
on small scales, it is not possible to derive its global prop- 
erties from individual measurements of the density in small 
regions: the characteristics of large-scale structure must be 
determined from a model of the fluctuating density field, and 
the values that one derives from it of statistical quantities, 
such as the mass variance. Similarly, the temporal variability 



of AGN is too large to ignore. We cannot hope to understand 
the radio luminosity function, or to test the AGN/cooling- 
fiow model, without a model of AGN variability. 

The normalized mechanical luminosities of a sample of 
AGN, L m = L m /Lx, should show scatter only because of 
temporal variability, and be free of differences in the scales 
of systems. Since the variability of AGN is known to be large 
and L m is an inherently positive variable, the natural variable 
to model is £ = ln(L m ), and the default assumption to make 
is that £(t) is a Gaussian random process. This assumption 
implies that individual measures of L m should have a log- 
normal distribution, so the time-sequence of L m is completely 
characterized by the power spectrum P(ui) of £{t). 

The enormous range of physical scales involved in AGN 
suggests that the power spectrum will be a power law P(uj) ~ 
a; - ' 3 . Several considerations suggest that /3 ~ 1, the index 
characteristic of flicker noise. First, on the short time-scales 
that are directly accessible, monitoring of AGN yields /3 ~ 1. 
Second, a comparison of the amplitudes measured in these 
studies with that inferred from the short duty cycle of quasars 
requires /3 ~ 1. Third, microquasars, whose power spectra 
can be probed over a camparatively wide range of frequen- 
cies, are found to have j3 — 1. Finally, diverse physical phe- 
nomena are found to be characterized by flicker-noise power 
spectra, so it would not be surprising if AGN conformed to 
this pattern. 

We estimate the spectral amplitude of the mechanical lu- 
minosity of AGN at the centres of cooling flows by assuming 
that energy radiated by the thermal plasma is on the aver- 
age replaced by the mechanical input of the central AGN, 
and that X-ray cavities probe the mechanical luminosity av- 
eraged over time-scales ~ 10 s yr. The amplitudes required 
are slightly larger than those that characterize the radio out- 
put of two microquasars, and of the order of those estimated 
from the duty cycle of quasars. 

We estimate the coefficients k in an assumed propor- 
tionality between the mechanical luminosity of an AGN and 
its radio power, either at 1.4 GHz or bolometric. We do this 
in two ways. In the first we directly fit a straight line to 
the correlation between radio and mechanical luminosities of 
sources that have radio- loud cavities. Alternatively, we fit our 
model of AGN variability to the radio data of all sources with 
cavities, whether radio-loud or ghost, under the assumption 
that (Lm) — Lx- The two methods yield compatible ranges 
of values of k. 

The available observations fail to constrain the model's 
parameters strongly, both because the currently available 
sample of sources is small and because our understanding of 
the physics of these sources is poor. Despite the significant 
advance that was made possible by the discovery of num- 
bers of X-ray cavities, measurements of the ratio of radio 
to mechanical luminosity are still extremely uncertain. Also, 
our analysis clearly shows the importance of determining the 
duty cycle of the radio sources in order to quantify their con- 
tribution in the energetics of cooling flows: in particular, it 
is crucial to estimate the amplitude of their power spectrum 
on time-scales >10Myr. We conclude that larger samples 
of AGN/cooling-flow pairs and better theoretical models of 
radio sources are needed to challenge the AGN/cooling-flow 
model. 

In a scenario in which AGN radio activity is strongly 
variable with time, the local luminosity function of radio 
galaxies should be interpreted as the distribution of the time 
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spent by a galaxy at a given radio luminosity rather than 
a distribution of sources with different intrinsic properties. 
This picture is supported by the observational finding that 
early-type galaxies hosting supermassive BHs with the same 
mass span several orders of magnitude in radio power. We 
showed that the FIRST-2dFGRS luminosity function of radio 
galaxies can be fitted by a log-normal distribution, consistent 
with the hypothesis that in all early-type galaxies, the cool- 
ing flow is balanced by a strongly variable radio source. Our 
model predicts a turnover at luminosities similar to the lowest 
currently probed. Any turnover in the luminosity function of 
AGN will be to a large extent masked by radio emission from 
nuclear starbursts. Nonetheless, the prediction of a turnover 
is interesting and will be tested by forthcoming radio surveys 
with instruments such as the Square Kilometer Array. 
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